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We study how deviations from spherical symmetry of a system, produced by angular momentum, and shear 
stress, modify the spherical collapse model, as linear density threshold for collapse of the non-relativistic com- 
ponent (6.) and its virial overdensity (Av), in Einstein-de Sitter and ACDM models. We modify the spherical 
collapse model to take account of the shear term and angular momentum term. We find that the non-spherical 
terms change the non-linear evolution of the system and that the collapse stops “naturally” at the virial ra- 
dius. Moreover, values of the linear overdensity parameter and of Ay are modified with respect to the standard 
spherical collapse model. 
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I. INTRODUCTION 


Current analyses of high quality cosmological data coming from Supernovae Type Ia (Riess et al. 1998; Perlmutter et al. 1999; 
Kowalski et al. 2008; Amanullah et al. 2010), CMB (Spergel et al. 2003; 2007), and clusters (Allen et al. 2002; Lima et al. 
2003; Eisenstein et al. 2005; Allen et al. 2008) are suggesting a cosmic expansion history involving some sort of dark energy and 
a flat spatial geometry in order to explain the recent accelerating expansion of the Universe. Among a number of possibilities 
to describe the dark energy (DE) component, the simplest one based on a cosmological constant A (see Padmanabhan 2003; 
Peebles & Ratra 2003; Lima & Braz 2004; Frieman et al. 2008; Li et al. 2011 for reviews), usually interpreted as the vacuum 
energy density (py) which acts on the Friedmann’s equations as a perfect fluid with negative pressure (py = —py). In the 
present cosmic concordance ACDM model, the overall cosmic fluid contains non-relativistic matter (baryons + cold dark matter, 
Qnr = 0.274) plus a vacuum energy density (Q4 = 0.726) that fits accurately the current observational data and thus it provides 
an excellent scenario to describe the present observed Universe (Komatsu et al. 2011; Allen et al. 2011; Del Popolo 2014). 
Nowadays, one of the most challenging problems in the so-called ACDM cosmology (see Del Popolo 2007, 2013, for a review) 
is to understand the role played by the different cosmic components during the non-linear regime of gravitational clustering and 
how the many possible physical effects contribute to determine the total mass of virialized halos (Del Popolo & Gambera 1996; 
Del Popolo 2002; Colafrancesco, Antonucio-Delogu, Del Popolo 1995; Del Popolo 2002) (galaxy and galaxy clusters)!. 

Most of the field of structure formation and galaxy formation concerns understanding the non-linear regime and collapsed 
objects in the density field. This can be done through N-body simulations or in some cases, using analytical models (Zel’dovich 
1965, 1970, Bisnovatyi-Kogan 2004, 2005, Bisnovatyi-Kogan, Tsupko 2015). 

A popular analytical approach to study the non-linear evolution of perturbations of dark matter (in the presence of a non- 
clustered dark energy (DE)) is the standard spherical collapse model (SSCM) proposed in the seminal paper of Gunn & Gott 
(1972) and extended in subsequent papers (Ryden & Gunn 1987; Gurevich & Zybin 1988a,b; White & Zaritsky 1992; Sikivie et 
al. 1997; Le Delliou & Henriksen 2003; Williams et al. 2004; Basilakos et al. 2010; Del Popolo & Kroupa 2009; Cardone et al. 
2011a,b; Del Popolo 2012a,b.). 
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' Note that the ACDM model suffers from other problems, like the small-scale problems (e.g., the cusp/core problem (Cardone & Del Popolo 2012; Del Popolo 
et al. 2013a; Del Popolo & Hiotelis 2014; Del Popolo, Le Delliou 2014, Del Popolo, Le Delliou 2017), and the missing satellite problem, the too-big-to-fail 
problem (Del Popolo & Gambera 1997; Del Popolo et al. 2014), and other problems like the cosmological constant problem (Weinberg 1989; Astashenok & 
Del Popolo 2012), and the cosmic coincidence problem. Moreover, is still debated the universal nature of dark matter density profiles forming in the ACDM 
cosmology (Navarro et al. 2010; Del Popolo 2010, 2011; Saburova, Del Popolo 2014; Del Popolo, Pace 2016). 


In the simplest case, termed “top-hat” model, one considers a uniform, spherical perturbation characterized by an overdensity 
ô = p(t)/p» — 1, where p, is the smooth background density of matter. In the quoted model, a proto-structure is considered as 
formed by concentric shells, expanding with the Hubble flow. The equation of motion for a shell in the perturbation is given by 


R=—-GM/R? (1) 
For a positively curved matter dominated universe the previous equation has the parametric form solution 


R œx (1-— cos(0)) 
t x (1-— sin(0)) (2) 


Starting from an initial comoving radius x;, each shell expands until 0 = m (t = tmax) then turn around and collapse at 0 = 27 
(t = 2tmax), formally reaching an infinite density?. 
At that time the density of the spherical region compared to the Einsten de Sitter (EdS) background is 
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On 
p/ Po = Te 7 5:55 (3) 


So the spherical perturbation starts to collapse when its density has reached 5.55 times of the background density. 
When the spherical region starts to collapse, the linear perturbation theory predicts that djin = 3(67)?/ 3/20 ~ 1.06. Assuming 
virialization to occur at t = 2tmax, the linear density contrast has at this point increased to 


2tmax 
bin = fe = 3/20 (6r ==)? ~ 1.686. (4) 


The model describes how a spherical symmetric overdensity decouples from the Hubble flow, slows down, turns around and 
collapse. In the last decade, the SSCM has been applied to study density perturbation evolution and structure formation in the 
presence of DE. However, when solving the density contrast (8) in the SSCM, the local shear (c) and the rotation (w) parameters 
are usually not taken into account. While the first assumption is correct, since for a sphere the shear tensor vanishes, the rotation 
term, or angular momentum is not negligible. A simple approach preserving spherical symmetry is to assume that the particles 
are described by a random distribution of angular momenta such that the mean angular momentum at any point in space is zero 
(White & Zaritsky 1992). Nevertheless, in any proper extension of the SSCM both effects need to be considered (Engineer et al. 
2000, hereafter E00; Del Popolo et al. 2013b,c,d) since shear induces contraction while vorticity induces expansion as expected 
from a centrifugal effect. 

In this paper, we study the net physical effect of shear and rotation in the framework of an extended spherical collapse model 
(ESCM). We restrict our analysis to the Einstein-de Sitter (EdS) and the flat ACDM background cosmologies. For the ACDM 
model we assume the following cosmological parameters: Qm = 0.274, Qa = 0.726 and h = 0.7. In particular, we discuss 
how the linear density threshold for collapsing non-relativistic component (6,) and its virial overdensity (Av) change. We recall 
that the change of these two parameters has a strong effect on the mass function and other fundamental cosmological quantities. 
As a general result, it is also found that the extra terms appearing in the ESCM are responsible for higher values of the linear 
overdensity parameter at galactic scales as compared to the case without shear and rotation. We also show that the non-spherical 
terms give rise to a collapse that “naturally” stops at the virial radius, differently from the SSCM, in which the collapse has a 
singular behavior predicting infinite density contrasts for all collapsed objects (see also Bisnovatyi-Kogan, Tsupko 2005, 2008). 
In real systems, the collapse to a point will never occurr in practice. Dissipative physics and the process of violent relaxation will 
eventually intervene and convert the kinetic energy of collapse into random motions leading the system to the virial equilibrium. 
The virial radius can be easily computed to be half the maximum radius reached by the system (ryir = 1/2r max). 

Even if the virialization argument is physically well motivated in real systems, in the SSCM no mechanism exists leading the 
system to virialization. Usually, one introduces by hand the assumption that in the collapse the shells constituting the system stop 
at a fixed radius (e.g., 1/2rmax). Taking into account angular momentum and shear in the SSCM, the system will not collapse 
to a point but the shells of the system will smoothly evolve towards fyir = 1/27 max. 


I. 6, AND Ay 


To begin with, let us now consider that the only clustering component in the cosmic medium is the cold dark matter. Following 
standard lines, the evolution of the overdensity ô in the SSCM, under the assumption that only DM can form clumps and that 


2 In this simple-minded model matter has no internal pressure, so there is nothing stopping the spherical blob to collapse to the infinite density. In real life 
collapse will, of course, stop before the infinite density is reached, giving rise to a “‘virialized” structure, when non-linear processes in the collapse phase 
convert kinetic energy into random motions. The final result will be a system which satisfies the virial theorem, and ryiy = 1/2rmax. 


DE is present as a background fluid (Fosalba et al. 1998; Ohta et al. 2003; Mota & van de Bruck 2004; Abramo et al. 2007), is 
given by a second order non-linear differential equation (Padmanabhan 1996; Ohta et al. 2003; Pace et al. 2010; Pace, Batista, 
Del Popolo,2014), namely 


5 . 4 §? 7 3 nan 


where the shear term o? = o;;0") and the rotation term w? = wj;j;w") are connected to the shear tensor, which is a symmetric 


traceless tensor, while the rotation is antisymmetric. 
Recalling that 5 = p/p — 1 = (a/R)? — 1 (ais the scale factor and R is the radius of the perturbation), and inserting it into 
Eq. 5, it is easy to check that the evolution equation for ô reduces to the SSCM (Fosalba et al. 1998; E00; Ohta et al. 2003) 


ËR rO GM ee 
TE = 4nGpR — 1/3(0° — w )R = -z — 1/80 —w*)R, (6) 


comparable with the usual expression for the SSCM with angular momentum (Peebles 1993; Nusser 2001; Zukin & Bertschinger 
2010): 


dR GM L? GM 4 
= i = QR? 7 
dt? R? M? R3 R2 = 25 ie 


where in the last expression we have used the momentum of inertia of a sphere, J = 2/5M R?. 

The previous argument shows that vorticity, w, is strictly connected to angular velocity, Q. In the simple case of a uniform 
rotation with angular velocity Q = Qez, we have that Q = w/2 (see also Chernin 1993, for a more complex and complete 
treatment of the interrelation of vorticity and angular momentum in galaxies). 

One assumption generally used when solving the SSCM equations for the density contrast ô (Eq. 5) is to neglect the shear, o, 
and the rotation w. While the first assumption is correct, since for a sphere the shear tensor vanishes, the rotation term, or angular 
momentum is not negligible. In fact, if we consider the ratio of the rotational term and the gravitational one in Eq. 7, we get 
a that for a spiral galaxy like the Milky Way, with L ~ 2.5 x 10’4gcm?/s (Ryden & Gunn 1987; Catelan & Theuns 1996) 
and the radius 15 kpc, is of the order of 0.4, showing, as well known, that the rotation is not negligible in the case of galaxy 
sized perturbations. The quoted ratio is larger for smaller size perturbations (of the dwarf galaxies size) and smaller for larger 
size perturbations (for clusters of galaxies the ratio is of the order of 10~°). The value of angular momentum, L, or similarly 
Q, can be obtained and added to the SCM as described in Del Popolo (2009) or as described previously, assigning an angular 
momentum x ,/GM(< r.)r, at turn-around (White & Zaritsky 1992; Sikivie et al. 1997; Nusser 2001). 

E00 studied the effect of the term o? — w? on the SSCM model just for the Einstein-de Sitter (EdS) model. 

In the present paper, we will study how the typical parameters of the SSCM (in Universes dominated by DE), namely ôe and 
Ay, are changed by a non-zero g and w terms. We also will show how the non-sphericity can “naturally” lead the system to 
virialization without the need to introduce it by hand, as usually is done in the SSCM. 

To this aim, we notice that Equation (5) can be written in terms of the scale factor a(t) as (Pace et al. 2010) 


n (3 F\ 94 6 3 Om E N 
(245 ô 3146 3 ab E(a) ot + 9) PH a) 7 w*)(1+6)=0. (8) 


The quantity Qm is the present-day value of the density parameter of the DM component while the quantity E (a) is defined by: 


Om 
Ba) = 4 -F + Qa, O) 


where Q; is the present-day value of the vacuum density parameter (at a = 1). 
In order to calculate the shear and vorticity terms in Eq. (8) it is convenient to define the dimensionless a-number as the ratio 
between the rotational and the gravitational term in Eq. (7): 


L? 


= MERG ` (10) 


3 As previously stressed, the non-trivial role of angular momentum in the SSCM has been pointed out in a noteworthy number of papers studying structure 
formation in DM dominated universes (see also Del Popolo, Ercan, Xia 2001; Del Popolo 2009; Zukin & Bertschinger 2010, Cupani et al. 2011). 
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As already stressed, the above quoted ratio, a, is of the order of 0.4, for a spiral galaxy like the Milky Way (L ~ 2.5 x 
104g cm?/s; R ~ 15 kpe (Ryden & Gunn 1987; Catelan & Theuns 1996a), larger for smaller size perturbations (dwarf 
galaxies size perturbations) and smaller for larger size perturbations (for galaxy clusters the ratio is of the order of 1076). 

Based on the above outlined argument for rotation one may calculate the same ratio between the gravitational and the extra 
term appearing in Eq. (8) thereby obtaining 


(0? —w*)Ho? = -> =O, (11) 


Note that this result is the same of that assumed by E00 (their Eq. 24) in the limit 6 >> 1. However, while the approximation 
given by Eq. (11) is good in the non-linear regime, as noticed by E00, it is insufficient to cover the larger range of density 
contrast (especially the quasi-linear regime) which is of interest to us. As done by E00, we expand (o? — w°) Hg ? in Taylor 
series and retaining the first two terms, we have 


3 AQm,0 
2 a 


(0? —w*)H)? = ô— AJ + B/8? (12) 

In particular, the values of A and B are calculated numerically (Spedicato et al. 2003) as in E00 by constraining the solution 
of the spherical collapse by means of the two-point correlation function or, as we did, comparing the threshold of collapse, ô< 
obtained by Sheth & Tormen (2002) which furnish the base to obtain the Sheth-Tormen mass function (see also Sheth & Tormen 
2001; Del Popolo & Gambera 2000; Hiotelis & Del Popolo 2006, 2013) with the 6, parameter which is obtained from Eq. (8). 
Note that the Sheth-Tormen mass function and the 6, from whom it is obtained using the excursion set theory and is in very good 
agreement with simulations (Yahagi et al. 2004; Del Popolo 2006). 


HI. RESULTS 


In this section we discuss some physical consequences of the ESCM discussed in this paper. In particular, we show how the 
non-spherical terms introduce physically virialization without the need to add it by hand (as previously stressed). 

Then we obtain the linear overdensity parameter 6. and the virial overdensity Ay. In Fig. la we show the result of the 
integration of Eq. (8), using Eq. (12) for w? — a”. We used the same units as in E00: y(x) = R/ryir and x = t/, with 
B = \/(8/3°)A/(GM)r3,.. The quoted equation is a non-linear equationn non integrable in terms of known functions. 
We can give a fit to Fig. 1a plot that can be usefull for readers, 


O pla? + p2x + p3 
at + qlr? + q2x? + qBx + q4 


y(x) 


(13) 


where p1 = 23.29, p2 = —14.74, q1 = —7.185,q2 = 41.35, q3 = —33.86, and q4 = 8.814 (coefficients with 95% confidence 
bounds) 

As expected on physical grounds, the function has a maximum and slowly decreases to unity for large values of xv. 
Tvir/ Rmax œ 0.6 and for t + oo, R > fyir. The value is slightly different from that obtained by E00 since in our model 
we have also the A term. 

In Fig. 1b we compare, similarly to E00, the non-linear density contrast in the modified SCM (dashed line) with that in the 
SSCM (solid line). Both quantities are plotted against the linearly extrapolated density contrast, d;. The plot shows that at the 
epoch corresponding to 6; ~ 1.686 the SSCM has, as known, a singular behavior while our model has a smooth behavior, with 
ô ~ 120 (value slightly different from that obtained by E00). When R = Ryax/2 Ozscm ~ 88. 

Since the deviations from spherical symmetry are smaller at early epochs and grow as the system evolves, the two curves 
should converge going towards smaller values of 6, as observed. However, since we used a Taylor expansion in 1/6, we cannot 
compare the two curves at 6 << 1. 

In the next plots, we show how 6, and Ay are changed by the non-spherical terms. 

In Fig. (2) (4 plots), we show the evolution of the linear overdensity parameter ô. (upper panels) and of the virial overdensity 
Ay (lower panels) for the same EdS and ACDM cosmologies. In the left panels, the analyses based on the ESCM are restricted 
to a halo of 10'! Mo /h since for galactic masses the effect will be enhanced, while on the right panels we consider also the 
effect of distinct masses. As before, we concentrate our analyses to three different mass scales: galactic (~ 101! Mo /h), groups 
(= 1018 Mo/h) and clusters (~ 10!° Mo /h). As expected from the analysis of Fig. 2, with the growth of the mass the effect 
of the extra term in the ESCM becomes negligible, and we recover the same values of the SSCM case. It is also worth to notice 
that the results for the ACDM model reduce to the ones of the EdS model for sufficiently high redshifts, since the influence of 
the cosmological constant becomes rapidly negligible. We will therefore concentrate only on the analysis of the left panels. For 
the different line colors and styles, we direct to the caption of the figure. 
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Figure 1: Left panel: plot of the scaled radius of the shell y(x) as a function of scaled time x (solid line). Right panel: non-linear density 
contrast in the SSCM (solid line) and in the extended ESCM (dashed line), plotted against the linearly extrapolated density contrast dz, 


As expected, the 6. for the ESCM is ~ 40% higher than for the SSCM case and it decreases towards high redshifts, since the 
effect of the extra term becomes smaller. For the EdS model, 6, decreases from a value of ~ 2.3 at z = 0 to ~ 2.1 at z = 10. As 
expected, the linear overdensity parameter for the ACDM model is smaller than the EdS one. This is understood by taking into 
account that, if we want to have the same number of structures now, we need to have a faster growth of structures to overcome 
the influence of the cosmological constant. This translates into a lower dc. 

In the lower panels we compare the behavior of Ay in the SSCM approach with the one predicted by the ESCM description. 
The red dashed (blue short dashed) curve show the standard and the extended results for an EdS model, while the green dotted 
curve represents a ACDM model. It is clear that the ESCM description affects also the virial overdensity parameter. In particular, 
we see that Ay is always constant in time for the EdS model. However, with the extra term its value increases reaching 
Ay * 185, about 4% higher than the standard result. The curve for the ACDM model approximates the EdS at high redshifts, 
as expected. Once again, higher masses are less affected by the ESCM correction term (lower right panel). 

The previous results have effects on the time and way galaxies form. As it should be expected, the extra term (a? — w?) 
slows down the collapse, and, as such, higher values for the initial perturbations are required in order to have a collapse 
at the same time of a spherical collapsing sphere. This has implication for mass function. 

The term co? — w suppresses, as expected, the high-mass tail of the mass function, since rare events are more sensitive 
to the background cosmology and to the collapse process. In general low-masses objects are not severely affected by 
rotation, but we should expect that for DE models different from ACDM we should observe a suppression in the number 
of objects already of several orders of magnitude for galactic masses. 

As we will describe in a future paper, we expect, in the case of the ACDM model that on small masses, the mass function 
is largely independent of the cosmological model, but it depends strongly on ô<. 

Since in the extended SCM this is higher, we should observe an increase in the number of objects at galactic scale at 
z = 0, to decrease at higher redshifts where the contribution of the non-linear term decreases. Moreover, we should 
observe a general decrement in the number of objects at high masses (up to M ~ 1014 Mo) to increase again to unity for 
masses of the order of 1015 Mo.This is explained by the fact that at such masses, the linear overdensity parameter are 
practically the same, therefore the mass function must not change. 

We conclude therefore that the term o? — w? has a strong impact on structure formation and that it is worth to 
investigate different parametrizations for the additional term. 


IV. CONCLUSIONS 


In this paper we have discussed how shear and rotation affect the standard spherical collapse model. The net effect of such 
quantities which is x (o? — w?) has been phenomenologically described by a power law on the density contrast depending on 
two parameters (A and B), fixed by comparing the threshold of collapse, 6, as discussed in Sheth & Tormen (2002), with the 6, 
value which is directly obtained from Eq. (8). We have focused our discussion on the influence of such an extra term on (a) the 
spherical collapse parameters 6, and Avy, and (b) the natural introduction of virialization in the SSCM. 
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Figure 2: Upper (lower) panels: evolution of ô. (Av) with respect to the redshift z for the EdS and the ACDM models. The left panels show the 
time evolution for both parameters at galactic scale (10'' Mo /h). The red curve represents the solution for the EdS model in the non-rotating 
case while the blue and the green curves represent the EdS and the ACDM model when rotation is included. On the right panels we compare 
the time evolution for three different masses, 10'', 101? and 10'° Mo/h. Different colors and line-styles correspond to different masses 
and different cosmological models: red dashed (orange dot-dashed) curve represents a halo of 10+! Mo/h in a EdS (ACDM) cosmological 
background, blue short-dashed (dot-short-dashed cyan) curve represents a halo of 10+? Mo /h for an EdS (ACDM) model, while the green 
dotted (magenta) curve stands for an object of 10'° Mg /h in an EdS (ACDM) model. 


The last point (b) was shown in Fig. (1), showing that the collapse does not go to a singularity, but reaches a maximum and 
goes softly down to the virialization radius (Fig. (1), left panel). Similarly the right panel of Fig. (1) compares the non-linear 
contrast in the SSCM and that in the modified one, ESCM, in terms of the linear density profile. The plot shows the different 
behavior of the SSCM and the ESCM. Concerning ôe (~ 40%), and Ay, as it should be expected, the extra term slows down 
the collapse, and, as such, higher values for the initial perturbations are required in order to have a collapse at the same time 
of a spherical collapsing sphere. It is also found that the extra term contribution is more important for galactic scales so that 
its contribution becomes negligible at high masses (galaxy clusters). In Fig. (2) we have numerically evaluated and compared 
the evolutionary behavior of both the ESCM and SSCM approaches. We have seen that both the linear and the non-linear virial 
overdensity in the extended spherical collapse model are enhanced with respect to the standard spherical case. Enhancements 
are more pronounced for ôe (~ 40%), while for Ay are only of the order of few percent. 

These results reinforce the importance of a more complete and rigorous treatment involving the effects of shear and rotation 
at the late stages of the collapsing halo history mainly for the galactic scales. 
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